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Abstract. We present a new few-parameter phenomenological model of light curves 
of eclipsing binaries and stars with ttansiting planets that is able to fit the observed light 
curves with the accuracy better than 1% of their amplitudes. The model can be used 
namely for appropriate descriptions of light curve shapes, classification, mid-eclipse 
time determination, and fine period analyses. 


1. Introduction 

Modern physical models of eclipsing binaries and transiting exoplanets (eclipsing sys¬ 
tems = ES) (e.g. Wilson & Devinney 1971; Wilson & Van Hamme 2014; Hadrava 2004; 
Prsa & Zwitter 2005; Prsa et al. 2011) are able to simulate their light curves with im¬ 
pressive fidelity. However, a solution of the reverse problem, to derive all important 
parameters of the model from observational data, is seldom unique. If we have obser¬ 
vational data of moderate or poor quality, we are often forced to diminish the number of 
free parameters and thus use simplified (and frequently physically inconsistent) models. 
A decision which parameters to fix or leave floating is usually nontrivial. 

Nevertheless, for several common practical tasks, such as: 1) determination of the 
mid-eclipse times, 2) light curve fitting, 3) description and classification of ES light 
curve shapes for the purpose of the current and future surveys like ASAS or GAIA, and 
4) fine period analysis, we do not need to know the detailed physics of the system. A 
good approximation of the shape of observed light curves is usually sufficient in such 
cases. 

We present here a few-parameter general phenomenological model of eclipsing 
system light curves in the form of the special, analytic, periodic functions that are able 
to fit an overwhelming majority of the curves with an accuracy better than 1%. Eor the 
sake of simplicity, we will limit our considerations only to systems with approximately 
circular orbits that represent more than 80% of observed eclipsing systems. 
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2 . The Model of Monochromatic Light Curves of Eclipsing Systems 

The model funetion of a monoehromatie light eurve (expressed in magnitudes) of eelips- 
ing systems F{(p, A), ean be assumed as the sum of three partieular funetions: 

F{<p, A) - F,{ip, A) + Fp(<f, A) + F,(<p, A), (1) 

where F^ deseribes the mutual eelipses of the eomponents, Fp models the proximity ef- 
feets, while Fc approximates the O’Connell effeet (irrespeetively of its physieal eause). 
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Figure 1. The Kepler light curve of an exoplanet transit in TrES-2b (Raetz et al. 
2014) is represented by points while the model ht is shown by a line. The differences 
between the observed magnitudes and the model, magnihed by a factor of 10, are 
shown in the lower part of the hgure. 


The normally occurring features of all eclipsing system light curves are two more 
or less symmetrical depressions caused by mutual eclipses of gravitationally bound 
stellar or planetary components. The model function of eclipses fe can be approximated 
by the sum of two symmetrical analytic functions of the phase (f. We assume that the 
primary eclipse is centered at the phase ^ = 0, while the secondary eclipse (if one is 
present) is centered at phase ip - 0.5, 
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ipk = [< 4 ^-{k- l)/2] - round [ip-{k- l)/2], 


where the summation is over the number of eclipses, - 1 or \ (the com¬ 

mon situation for exoplanet transits), (pk are the auxiliary phases defined in the interval 
(-0.5,0.5), Ak are the central depths of eclipses in mag, is the parameter express¬ 
ing half-widths of eclipses, Ck,K are the correcting parameters, /]. parameterize the 
kurtosis of individual eclipses. 
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Figure 2. V light curve of an unusual eclipsing binary V477 Lyr consisting of a 
very hot nucleus of a planetary nebula and solar type star. Only eight free parameters 
is needed for the ht of the light curve strongly affected by proximity effects. V data 
taken from Pollacco & Bell (1994). 


The fit of a transit of the exoplanet TrES2b (tig - 1) observed by Kepler (see 
Fig. 1) (Raetz et al. 2014) utilizes five parameters: A\,C\,K,D\,ri. However, the full 
deseription of the model monoehromatie light eurves with two eelipses (tig = 2) would 
generally require nine free parameters: Ai 2 , Ci^ 2 > ^i, 2 > A, 2 > and K. Nevertheless, the 
most common situations would permit a reduction to 5 parameters, Ai 2 ,C,D,r, be¬ 
cause it is usually acceptable to assume that Ci C 2 , K - 0, Di - D 2 , Ei - r 2 . A 
detailed discussion of the properties of the eclipse model function Fg and its possible 
simplifications are given in Mikulasek (2015). 

Light variations of eclipsing binaries caused by eclipses are usually modified by 
proximity and O’Connell effects that may be modelled by a sum of cosine and sine 
harmonic polynomials: 


Hp+rie Hc+rip+ne 

Fp{(p) = ^ AkCos[27r{k-ns)(p]-, Fc(^) = ^ A*; sin(27r^), (3) 

^=ne+l k=np+ne+l 


where k = + \,nf. + 2,... ,{ni, + Up + n^), Up is the number of terms in Fp(^i): 

Hp = 0,2,3,..., Wg = 0 if the O’Connell asymmetry is not present^ else \ (for 
more information see Mikulasek 2015). 

The uncommon V light curve of eclipsing binary V477 Lyr (see Fig. 2) is deter¬ 
mined by 8 parameters: Ai 2 , 3 , 4 , 5 , D, C, E, Up = 3, = 0. 


'If p > ^ then hi, = 0. 
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3. Multicolour Light Curves of Exoplanet Transits and Eclipsing Binaries 

The parameters of the above defined model funetions, espeeially the amplitudes Ak{A), 
and parameters Ck{A), D{A), and r{A) are generally wavelength dependent. It follows 



Figure 3. Light curves of the 10-colour photometry of the exoplanet transit of 
HD 209458b. Effective wavelengths of individual colours are in nm. The data were 
taken from Knutson et al. (2007) 


both from the theory and experienee that most of the wavelength dependeneies ean be 
well approximated by the low order polynomials of the quantity A; A = To/Teff - 1, 
where deft is the effeetive wavelength of eorresponding passband, Tq is an arbitrarily 
seleeted eentral wavelength of the data set. Then: 

Sak gel gd gy 

A, = 2]a,yA^-\ Ci = Y,cijA^-\ D = Y,djAj-\ (4) 

y=i ;=1 i=i i=i 

where gak,gci,gd,gy, I - 1,2 or Z = 1, are the numbers of degrees of freedom of 

the eorresponding parameters of the model. The standard set of the monoehromatie 
light eurve model parameters of eelipsing binaries (see See. 2): {A^, C/, D, T) ean be 
eonsidered as the speeial ease of the multieolour deeomposition Eq. 4: {aki,cii,di,yi}. 

All needed parameters and their uneertainties ean be determined using the standard 
non-linear weighted least square teehnique (see in Mikulasek 2015). 
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Figure 4. BVRI light curves of a ‘cool Algol’ AV Del, consisting of a main se¬ 
quence F type primary and the K subgiant filling its Roche lobe, can be approximated 
by our model function with 9 free parameters. The data were taken from Mader et 
al. (2005). 


4. Conclusions 

The outlined phenomenologieal modeling of eelipse-binary and planetary-transit light 
eurves has an advantage of high effieieney and simplieity and thus ean be used for a 
wide range of applieations. 
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